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Quantum computers could perform certain tasks which no classical computer 
can perform in acceptable times. Josephson junction circuits can serve as 
building blocks of quantum computers. We discuss and compare two designs, 
which employ charge or magnetic flux degrees of freedom to process quan- 
tum information. In both cases, elementary single- qubit and two-qubit logic 
gates can be performed by voltage or flux pulses. The coherence time is long 
enough to allow a series of such operations. We also discuss the read-out, 
i.e. a quantum measurement process. In the charge case it is accomplished 
by coupling a single- electron transistor to the qubit. 

PACS numbers: 85.25.Cp, 85.25.Hv, 03.67.Lx 



1. INTRODUCTION 

Quantum computers, if available, could perform certain operations in 
a massive parallel way, which would lead to an exponential speed-up as 
compared to the performance of classical computers. Furthermore, medium- 
size quantum information processing devices could be used in a number of 
applications such as, for instance, quantum communication and quantum 
cryptography. Several physical realizations of quantum bits, i.e. the elemen- 
tary building blocks of quantum computers, have been proposed. The best 
studied ones are trapped ions, NMR in the liquid state and cavity QED. Here 
we discuss solid-state nano-electronic realizations. They appear promising 
since they can be scaled up to large numbers of qubits and are most easily 
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embedded in electronic circuits. 

Josephson junction circuits are particularly suitable for quantum in- 
formation processing, since they combine the intrinsic coherence of the su- 
perconducting state and the possibility to control the circuit dynamics by 
voltage and magnetic flux pulses. Their fabrication and manipulation are 
possible by present-day technologies. The dynamical variables in these cir- 
cuits are the charges on the islands and the phases of the superconducting 
order parameter. Both are canonically conjugated, and Heisenberg's uncer- 
tainty principle holds for them. Depending on the ratio of two characteristic 
energies - the typical charging energy, which favors well-defined charges, 
and the Josephson energy, which favors the phase degree of freedom - the 
charge or the flux can have a well-defined value, while the other fluctuates 
strongly. Either chargeSi'i or phaseBi degrees of freedom can be used to 
store and process quantum information. Here we describe charge and phase 
(flux) nano-electronic quantum bits and possible designs of the circuits in 
both cases. 

The basic elements of a quantum computer are the qubits, i.e. two-state 
quantum systems which can be manipulated, separately for each qubit, by 
the control of the Hamiltonian. This allows performing single-bit logic gates. 
Additionally, one needs to be able to couple qubits in a controlled way to 
perform two-bit logic operations. The whole system should stay coherent 
for a sufficiently long time, that is, the coupling to the environment needs 
to be weak. Finally, to read-out the information about the quantum state of 
the system a quantum measurement should be performed. Below we discuss 
these steps for the relevant designs. 

2. JOSEPHSON-JUNCTION QUANTUM BITS 

2.1. Charge qubit 

The simplest design of the Josephson qubit using the charge degree of 
freedom is presented in Fig. [l]a. It consists of a superconducting electron 
box with a low-capacitance Josephson junction, with capacitance Cj and 
Josephson energy Ej, biased by a voltage source through a gate capacitor 
C g . If the superconducting gap is large enough then at low temperatures 
odd-parity states are forbidden,! and the charge on the superconducting 
island is a multiple of the Cooper-pair charge 2ne, measured relative to the 
neutral state. The Hamiltonian of the system is then given by 
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Fig. 1. a. The simplest design of a charge qubit. b. A charge qubit with 
the flux-controlled Josephson coupling, c. The rf-SQUID, the simplest flux 
qubit. d. A flux qubit with controlled barrier height between the wells. 



where the phase difference across the junction, ip, is canonically conjugated 
to n. The charging energy, with scale Eq = e 2 /2(C g + Cj), is chosen to 
dominate over the Josephson term Ej. In equilibrium at low temperature, 
k^T <C Eq, the system is in the ground state, which - away from certain 
degeneracy points - is approximately the charge state with the minimal 
charging energy. Only near the voltages Vd eg = (2n + l)e/C g , where the 
states with n and n + 1 Cooper pairs on the island are degenerate, the 
weak Josephson coupling mixes the charge states strongly. Biased near these 
voltages the system reduces to a two-state problem, with the Hamiltonian 
in the basis of charge states \n) and \n + 1) 

U= £ -a z + ja x . (2) 

Here e(y x ) = 2e^j-(V x — V^ cg ) denotes the difference in charging energy 
between two relevant charge states, the tunneling amplitude between the 
states is A = — Ej, and the capacitance of the qubit in the circuit is C q ^ = 

The first term in the Hamiltonian (0) can be controlled through the 
gate voltage. The Josephson coupling in the second term can also be con- 
trolled if the junction is replaced by the dc-SQUID threaded by a mag- 
netic flux <3? x , as shown in Fig. [l]b. Then the effective Josephson energy is 
£j($ x ) = 2£^cos(27r$ x /$o)- With the control over these two parameters 
one can perform any single-bit logic operation. In the idle state between the 
operations one keeps the qubit at the degeneracy point V x = Vd eg and chooses 
= c ^o/2- Hence the Hamiltonian vanishes, TC = 0, and the qubit's state 
does not evolve in time. To perform an operation, one can switch the flux 
to a different value for a finite time r. The resulting change of the quantum 
state of the qubit is described by the unitary operator ex.p(ia x Ejt/2). A 
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voltage pulse results in another elementary operation, exp(— ia z e(V x )t/2). 
With a series of (no more than three) such operations with proper time 
spans any 2x2 unitary operation, that is any single-qubit logic gate can be 
performed. The typical times of the operations are of order h/A or h/e. 

The manipulations of the state of the qubit are similar to the manipula- 
tions used in NMR experiments, and various techniques familiar from NMR 
applications can be employed. For instance, instead of rectangular voltage 
pulses one can use resonant ac-pulses to induce coherent transitions between 
the qubit's states. Then the typical operation time is determined by the 
amplitude of the ac-pulse and can be optimized. It should be slow enough 
to make the control easy and fast enough to maximize the number of logic 
operations performed within the phase coherence time. 



A controllable two-state quantum system can be realized also in the 
opposite limit of dominating Josephson coupling. The simplest design is 
the rf-SQUID (see Fig. [l|c) that is a superconducting loop interrupted by a 
Josephson junction. The phase difference across the junction, 2tt&/&q, is 
controlled by the flux $ in the loop, which fluctuates around the externally 
applied value <3? x . With the Josephson, charging and magnetic contributions 
taken into account, the Hamiltonian of the system reads 



Here L is the self-inductance of the loop and Cj the capacitance of the 
junction. The charge Q = —(i/h)d/d<& on the leads is canonically conjugated 
to the flux <E>. If the self-inductance is large (0l = Ej/(^1/4tt 2 L) is slightly 
larger than 1) and the externally applied flux is close to &o/2, the two first 
terms in the Hamiltonian form a double- well potential near $ = Qq/2. At 
low temperatures only the lowest states in the two wells contribute to the 
physics of the system. The reduced Hamiltonian of this two-state system is 
again given by Eq.(|), where now £(<3? x ) = 4tt^6(Pl — l)-Ej($x/$o - 1/2) is 
the asymmetry of the double well potential, and A is the tunneling amplitude 
between the wells. The latter can be controlled through the height of the 
barrier, which is determined by Ej. This Josephson energy can be controlled, 
in turn, if the junction can be replaced by the dc-SQUID, as shown in Fig. [l|d. 

With two external parameters governing the Hamiltonian, elementary z- 
and x-rotations can be performed, as we have seen in the previous subsection. 
They can be driven either by switching the external flux for a finite time or 
by resonant pulses. 

The rf-SQUID described above was discussed in connection with the 
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'macroscopic quantum coherence', that is coherent oscillations of a quan- 
tum system between two macroscopically different states. BH However, the 
requirements of sufficiently large self-inductance and Josephson energy of 
the junction make the rf-SQUID very susceptible to external noise, and the 
experiments with the rf-SQUID were not successful so far. To overcome this 
difficulty Mooij et al.El suggested to use a smaller superconducting loop with 
three junctions (one with controllable critical current). Then the double- well 
potential is formed by Josephson terms for the junctions, and lower critical 
currents can be used. As a result the system should stay coherent for a 
longer time. We shall discuss dephasing effects in Section |j. 



3. COUPLING OF THE QUBITS 

To perform a quantum computation with a register of qubits two-bit 
logic operations are necessary. For such an operation, which is a unitary 
transformation of the quantum state of two qubits, the couplings between 
the qubits should be controlled individually for each pair. In this section we 
discuss realization of such interactions for charge and flux qubits. 

3.1. Controlled coupling of charge qubits 

One possibility to couple charge qubits is to join them in parallel to a 
common inductor, as pictured in Fig. ^a. Then their dynamics is coupled to 
the oscillations in the LC-circuit, which is formed by the inductor and the 
capacitances of all N qubits in parallel. If the frequency u>lc = W LNC q ^ 
of these oscillations is higher than typical qubit frequencies, the oscillatory 
degrees of freedom are not excited by the qubit manipulations, but still they 
provide the effective coupling between the qubitsBH 

To clarify the physics of the coupling and to estimate its magnitude, 
we provide a simple derivation. The current through the inductor, /, is 
given by the contributions of all qubits, I = 5^/*. In terms of the flux 
$ through the inductor the Hamiltonian of the oscillations is expressed as 
H osc = g-/$ + ^or 

($-LI) 2 Q 2 LP 
Hosc ~ 2L + 2NC qh 2 ' {) 

Here Q = —(i/h)d/d<& is the charge canonically conjugated to <£. If the 
time-evolution of the qubits is slow compared to the oscillations, we can 
use an adiabatic approximation and consider I as constant, and the os- 
cillator remains in its ground state. The energy of the ground state of 
the first two terms in (||), %ujlc/2, does not depend on /, while the last 
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Fig. 2. Coupling of qubits via an oscillator circuit: a. for charge qubits and 
b. for flux qubits. 

term provides the current-dependent correction. This correction describes 
the effective coupling between the qubits, — L(J2i h) 2 /^- In the basis of 
the qubits' charge states individual current operators can be expressed as 
P = (Cqb/Cj)(eEj /h)(Ty, leading to the interaction term 
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This coupling can be controlled through external fluxes which bias the 
SQUIDs of the qubits. Between the operations we keep all $^ = 3*0/2, 
so that Ej = and the coupling is off. During single-bit operations we 
switch on Ej for at most one qubit, so that the interaction term is still zero. 
When a two-qubit operation is needed, the fluxes are changed for two qubits, 
providing the interaction between them for a finite period. This interaction 
leads to a non-trivial unitary transformation of the quantum state of two 
chosen quantum bits, i.e. to a two-bit logic gate. 

This design allows performing both single- and two-bit operations by 
controlling the same external parameters (fluxes and voltages) as for indi- 
vidual qubits. 

3.2. Coupling of flux qubits 

One possibility to couple flux qubits is to use a flux transformeJl which 
provides an inductive coupling between the qubits. Each qubit has two 
loops, those of the rf- and dc-SQUIDs in the simplest rf-SQUID design. Any 
loop of one qubit can be coupled to any loop of the other, giving rise to 
<7*<t|, a l z a J x or a l x a J x coupling terms. To turn off this coupling completely, 
one would need to have an ideal switch in the flux transformer. This switch 
is to be controlled by high-frequency pulses, and the related external circuit 
can lead to decoherence effects. An alternative is to keep the interaction on 
constantly and use ac driving pulses to induce coherent transitions between 
the levels of the two-qubit system (cf. Refs. [HH). A disadvantage of this 
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approach is that permanent couplings lead to unwanted 2-qubit correlations 
between the logic operations. 

Here we discuss another design of a many-qubit circuit, which allows to 
control interactions via flux sources of individual bits. The circuit, shown 
in Fig. |||b, is similar to the register of charge qubits discussed in the previ- 
ous subsection. It includes an LC-circuit formed by a loop, with the self- 
inductance L osc , interrupted by a small capacitor C osc . The loop is coupled 
inductively to the set of flux qubits and mediates interaction between them. 
Following the derivation in Refs. |], ^ one can find the effective coupling by 
integrating out the oscillations which are faster than the qubit's dynamics. 
A simple way to obtain the coupling is to notice that in the low-capacitance 
limit, Cqsc - > (almost decoupled qubits) the effect of the qubits is to estab- 
lish a voltage drop across the inductor L osc , given by V = J2i M&i/L. Here 
for each qubit $j is the flux in its loop, L is the self-inductance of the loop 
and M is the mutual inductance with the loop of the LC-circuit. Then the 
Hamiltonian for the oscillator mode is Ti. osc = $ 2 /2L osc + Q 2 /2C osc — VQ, 
with Q being the charge conjugated to the flux <3? through the LC-circuit. 
By similar arguments as in the previous subsection we find that the oscil- 
lator provides the inter-qubit interaction term — C osc V 2 /2. In limit of the 
weak coupling to the LC-circuit, the time derivatives of the qubits' fluxes are 
given by <&i = ^[Hi,^%] = 5$>iAi(Ty/h where 5&i is the separation between 
two minima of the potential of the rf-SQUID. Hence, the interaction is given 



This interaction can be controlled through tunneling amplitudes Aj of in- 
dividual qubits. To turn off the interaction, one should put Aj to zero. 
The amplitudes can be exponentially suppressed by increasing the barrier 
heights. (Note that in this case also the fluctuations of Aj are exponen- 
tially suppressed, unlike for the charge qubits, where stronger fluctuations 
are present.) The absence of the interaction between the operations helps 
increasing the accuracy of the manipulations. 



Due to the unavoidable coupling to environmental degrees of freedom 
the quantum state of the qubits gets entangled with the environment, which 
leads to dephasing effects. Charge qubits are sensitive to the electromagnetic 
fluctuations in the external circuit and in the substrate and to background 
charge fluctuations. Flux qubits are insensitive to the latter. Various sources 
of dissipation for flux qubits are discussed in Ref. || Here we estimate the 
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4. ENVIRONMENT AND DEPHASING 
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effect of fluctuations in the external circuit. 

The control over the voltage or flux requires coupling the quantum sys- 
tem to a dissipative external circuit, which introduces fluctuation effects. 
The voltage fluctuations in charge qubits or the flux fluctuations in the 
main loop of flux qubits are coupled to cr 2 -degree of freedom. The external 
circuit can be parameterized by an effective impedance of the voltage (flux) 
source, which in turn can be modeled by an oscillator bath in the spirit of 
the Caldeira-Leggett model. The behavior of a two-state system coupled 
to an oscillator bath has been discussed in the context of the spin-boson 
modeller O The conclusion is that the decay of the off-diagonal elements 
of the density matrix of the qubit (dephasing) and the relaxation of the di- 
agonal elements to their equilibrium values are described by the two time 
scales 

h (1 AE 2 2k B T 2 X- 1 

h ( , ae , , r 1 

W = 7 ^ { coth ^f sm > (8) 

respectively. Here h/AE gives the typical operation time (cf. Section 2) in 
terms of the level spacing AE = Ve 2 + A 2 of the qubit, T is the temperature, 
and tan 77 = A/e. 

The dimensionless parameter 7 describes the strength of the dissipation. 
For charge qubits it is determined by the resistance Ry of the voltage source 
in units of the large quantum resistance Rk = h/e 2 ~ 25.8kf2, 

A small gate capacitance C g ~ C q b <C Cj provides a weak coupling between 
qubit and environment. 

For flux qubits, 7 is fixed by the impedance Rj of the current source in 
the input loop, which provides the flux bias, 

v-^K*) ■ <10) 

Here M is the mutual inductance of the input loop and the qubit 's loop. The 
numerical prefactor in Eq. (^o|) is a = 7t/(6(/3l — 1)) for the rf-SQUID and is 
also of the order of unity for the design of Ref. ||. The dephasing is slow for 
small loops and junctions with low critical currents. Indeed, the argument 
in the bracket in ( |io| ) is proportional to $q/47t 2 L£'j. While this quantity 
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should be slightly smaller than one for the rf-SQUID, it can be much larger 
for the design of Mooij et al.,i corresponding to slower dephasing. 

The qubit is equivalent to a spin- 1/2 particle in the external magnetic 
field ez + Ax, while the environment produces an additional fluctuating field 
in z-direction. The component of the fluctuating field orthogonal to the 
external field (with the magnitude proportional to sin 77) induces transitions 
between the eigenlevels, while the longitudinal component (cc cos rj) leads 
to random fluctuations of eigenenergies. Only the former process leads to 
relaxation (|8|), while both contribute to the dephasing ([?]). This observation 
explains the ^-dependence of the decoherence rates (0), (|8|). The relaxation 
rate is slow for low-impedance voltage sources and high-impedance current 
sources. Weak coupling to the external circuit, i.e., small C g or small M, 
helps further maintaining the coherence. 

The effect of fluctuations of the a x -term in the Hamiltonian (in the flux 
circuit of the dc-SQUID-loop which controls the Josephson coupling in both 
charge and flux qubits) can be described in a similar wayJl Because of the 
different direction of the fluctuating effective 'magnetic field', sin 77 and cos 77 
should be interchanged in those terms. The effect of these terms is relatively 
weak for the operation regimes discussed in Refs. [2] and ||. 

5. READ-OUT: QUANTUM MEASUREMENT 

To complete the quantum computation, one needs to read out the infor- 
mation about the quantum state of the qubits, that is to perform a quantum 
measurement. For charge qubits this can be accomplished by coupling the 
qubit capacitively to a single-electron transistorH3 as shown in Fig. ||. During 
the computation the SET is kept in the off-state Vt T = 0, with no dissipative 
current and no additional decoherence. The only effect of the transistor is 
a renormalization of the capacitances in the qubit's circuit. To perform the 
measurement, the transport voltage is switched to a sufficiently high value 
and the dissipative current starts to flow. The value of the current in the 
circuit of the transistor is very sensitive to the charge on the island of the 
qubit. Monitoring the current, one can extract from the data information 
about the state of the qubit. 

To study this quantum measurement process, we analyzed the time 
evolution of the density matrix of the coupled system of the qubit and the 
SET0! The system is characterized by three energy scales: the typical 
Coulomb energy of the transistor, -E se t, the charging energy of the qubit, 
e(V x ), and the Coulomb interaction between the charges of the qubit and the 
middle island of the SET, E mt . We choose E set to be the largest energy scale, 
E sct » f » Ei nt . The Josephson coupling Ej($ x ) is switched to a small 
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Fig. 3. Read-out is accomplished by coupling a single-electron transistor 
capacitively to the qubit. 



(or zero) value before the measurement. The transport voltage should be 
large enough, of order E se t/e, to overcome the Coulomb energy gap between 
two charge states of the SET. The value of this gap is slightly different for 
different states of the qubit, which leads to different tunneling rates r ± 5T 
through the transistor, with T of order 2TraE set and ST of order 2iraEi nt . 
Here 2ira = Rk/^ 2 Rt <C 1 is the dimensionless tunneling conductance of 
the junctions of the transistor, of the same order in both junctions. For 
definiteness we consider the limit E[ nt <C T <C AE, where essential features 
of the evolution can be seen, although the derivation can be performed in a 
wider parameter range. 

To describe the system we single out the most important degrees of 
freedom: the qubit's degree of freedom, the charge Ne on the middle island 
of the SET and the charge me which passed through the transistor after 
the transport voltage is turned on, that is the time integral of the current. 
The microscopic degrees of freedom can be traced out, and a closed set of 
equations can bederivedBEl for the elements of the density matrix, p*- (N, m) 
which are diagonal in N and m. 

The evolution of the system is characterized by three time scales. In 
the first stage the random tunneling processes in the SET lead to the loss of 
the phase coherence of the qubit. The off-diagonal elements of the reduced 
2x2 density matrix of the qubit, J2n mPj(N,rn), decay to zero on a short 
time scale of order t v « Y / Ef nt . 

At a later stage the information about the qubit's state can be deduced 
from the current. The dynamics of the current in the SET can be described 
by the probability distribution of the number of electrons m which have 
passed through the transistor: P(m,t) = J2i n P\(N, m)- At the start of the 
measurement t = no electrons have tunneled, and P(m,t) = 5 m fi. Then 
the current begins to flow, and the peak shifts in the direction of positive 
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Fig. 4. The evolution of the probability distribution p(m/t,t) for the values 
of the average current in the SET. The initial state of the qubit is a = ^/l/4, 
b = y/3/4. a. t < t mcas : a very broad peak due to shot noise, b. i mea s < 
t < t m i x : two peaks separate, c— e. t ~ i m i x and t > f m i x : a plateau grows 
between the two peaks, erasing the information about the initial state of the 
qubit. 

m and widens due to shot noise effects. Since the two charge states of the 
qubit correspond to different net currents in the SET, after some time the 
peak splits into two. The separation of the peak centers grows linearly with 
time, 25Tt, while their width grows as yTt. Therefore they separate after 
time t m eas ~ T/5T 2 . The weights of the peaks after the separation are given 
by |a| 2 and \b\ 2 for the initial state of the qubit a\0) + b\l), and a good 
quantum measurement is realized if m is measured after t meas . As expected, 
the measurement time is longer than the dephasing time. 

This ideal picture gets more complicated at even longer times. At fi- 
nite, or fluctuating, Ej the eigenstates of the qubit are different for different 
charges Ne of the SET, and the measurement induces transitions between 
them. After a long mixing time t mix AE 4 /(E 2 nt E%T) the diagonal el- 
ements of the density matrix of the qubit become equal to 1/2, and the 
information about their initial values is lost completely. At the same time 
P{m, t) develops a plateau between the two peaks. The weights of the peaks 
decay exponentially. At longer times t 3> i m ix the plateau transforms into a 
peak around m = Tt. This distribution contains no information about the 
initial state of the qubit anymore. 

The time evolution of P(m, t) is depicted in Fig. |I[ For convenience the 
probability distribution p(m/t, t) of possible values of the quantity m/t, the 
current in the SET averaged over time t, is shown. It is related to P(m, t) 
by p(I, t) = t P(It, t) where the prefactor t ensures the normalization of the 
distribution. If the current is averaged over short times t < i mea s) Fig. ||a, 
the shot noise does not allow to distinguish between close values of the 
current corresponding to different qubit states (very broad peak around er). 
At longer times i mea s < t < i m i x , Fig. ||b, two peaks are formed around 
eT± = e(T ± ST), and the measurement can be performed. At t ~ i m i x a 
plateau starts to grow between the peaks, and at the very long times t 3> t m j x 
the plateau takes over and transforms into a narrow peak around eT. 



Yu. Makhlin, G. Schon, and A. Shnirman 



6. CONCLUSION 

We have discussed two possible designs of nano-electronic quantum bits, 
based on Josephson junction circuits in the charge and flux regime. In both 
cases the quantum dynamics of the qubits can be controlled through voltages 
and fluxes (currents) . We suggest a way to couple flux qubits which is a dual 
analog of the approach suggested earlier for charge qubits. Compared to the 
proposal of Mooij et al.B it has the advantage that it allows to control the 
inter-qubit interaction without introducing new links between the qubits 
and the external circuit. Apart from that we estimated the dephasing times 
due to coupling to electromagnetic fluctuations in the external circuit. This 
contribution has not been discussed earlier for flux qubits. The quantum 
measurement of the state was discussed in detail for charge qubits. We 
provided an analysis of the long-time evolution of the probability distribution 
of current values in the single-electron transistor and of the qubit's density 
matrix. The design of a quantum measurement circuit for flux qubits requires 
further investigation .0 

Earlier experiments with the superconducting box have demonstrated 
the quantum nature of the two-state system and superpositions of different 
charge states. Bill Recently Nakamura et al.0 observed time-resolved co- 
herent oscillations of the quantum state of a charge qubit. Their experiment 
is the first observation of the macroscopic quantum coherence (MQC) effect 
and realizes a single-bit operation (corresponding to an x-rotation). The 
observed phase coherence time of several nanoseconds was shorter than the 
value ~ 100 ns, predicted by (0), (|), (§) 

because of the use of a simple (but 
efficient) measuring procedure. Further development of the measurement 
apparatus (cf. Ref. |l7|) can render the coherence time longer. In flux sys- 
tems, the attempts to observe MQC with the rf-SQUID were not successful, 
but recent developments} should make it possible in the near future. 
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